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(57) ABSTRACT 

A system for receiving, analyzing and communicating results 
of sensing chemical and/or physical parameter values, using 
wireless transmission of the data. Presence or absence of one 
or more of a group of selected chemicals in a gas or vapor is 
determined, using suitably functionalized carbon nanostruc- 
tures that are exposed to the gas. One or more physical param- 
eter values, such as temperature, vapor pressure, relative 
humidity and distance from a reference location, are also 
sensed for the gas, using nanostructures and/or microstruc- 
tures. All parameter values are transmitted wirelessly to a data 
processing site or to a control site, using an interleaving 
pattern for data received from different sensor groups, using 
I.E.E.E. 802. 1 1 or 802. 15 protocol, for example. Methods for 
estimating chemical concentration are discussed. 

10 Claims, 13 Drawing Sheets 
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NANOSTRUCTURE SENSING AND 
TRANSMISSION OF GAS DATA 

ORIGIN OF THE INVENTION 

Work related to this invention was performed under a joint 
research agreement between the National Aeronautics and 
Space administration and the University of California. 

FIELD OF THE INVENTION 

This invention relates to wireless transmission of data pro- 
vided by nanostructure-based chemical and physical sensors. 

BACKGROUND OF THE INVENTION 

Chemical sensors have been developed for decades to 
detect various concentration levels of gases and vapors for 
deployment in a wide range of applications in industry, space 
mission, environment monitoring, medical, military, and oth- 
ers. The detection usually centers on change in a particular 
property or status of the sensing material (such as thermal, 
electrical, optical, mechanical, etc) upon exposure to the spe- 
cies of interest. The sensing material may be any of several 
elements from the periodic table, plus inorganic, semicon- 
ducting and organic compounds, in bulk or thin film form. 
One of the most investigated classes of chemical sensors in 
the past is the high-temperature metal oxide sensors due to 
these sensors’ high sensitivity at low ppm to ppb concentra- 
tion levels, with tin oxide thin films as an example. Polymer 
sensors have been studied in recent years because they can be 
operated at room temperature with low power consumption 
and are easily fabricated. Although commercial sensors based 
on these materials are available, continued research is in 
progress with sensing technologies using new sensing mate- 
rials and new transducer platforms. New sensing technolo- 
gies, such as nanotechnology -based sensors, are being devel- 
oped to overcome the large power consumption and poor 
selectivity of metal oxides sensors, and to improve the poor 
sensitivity and narrow detection spectrum of polymer sen- 
sors. 

Typical figures of merit expected from a chemical sensor 
include sensitivity (even down to a few molecules, selectivity, 
low power consumption, rapid response time and rapid sensor 
recovery time. Sensors based on the emerging nanotechnol- 
ogy may provide improved performance on all of the above 
aspects compared to current micro and macro sensors. Nano- 
materials exhibit small size, light weight, very high surface- 
to -volume ratio, and increased chemical reactivity compared 
to bulk materials; all these properties are ideal for developing 
extremely sensitive detectors. The potential of nanomaterials 
for detection and for protection and remediation has been 
outlined in a recent report by the National Science and Tech- 
nology Council (“Nanotechnology Innovation for Chemical, 
Biological, Radiological and Explosive Detection and Pro- 
tection,” November 2002). 

One promising nanomaterial is the carbon nanotube 
(“CNT”), which exhibits extraordinary mechanical, electrical 
and optical properties. These interesting properties have 
prompted wide range investigations for applications in nano- 
electronics, high strength composites, field emitting devices, 
catalysts, etc. A single-walled carbon nanotube (“SWCNT”) 
has all the atoms on the surface and therefore would be 
exposed maximally to the environment, allowing a change in 
its properties sensitively. The first demonstration of an 
S WCNT-based sensor was for a chemical field effect transis- 
tor (“CHEMFET”), where a single semiconducting SWCNT 


2 

is used as the channel material and the conductivity was 
shown to change upon exposure to N0 2 and NH 3 (J. Kong et 
al, “Nanotube Molecular Wires as Chemical Sensors,” Sci- 
ence, vol. 287, pp 622-625, January 2000). The potential for 
5 using CNTs in chemical sensors was noted in the Kong et al 
article. However, it still is a challenge to make practical sen- 
sors, due to difficulty in fabrication complexity, low sensor 
yield and poor reproducibility. 

A different configuration of carbon nanotube-based chemi- 
io cal sensors with much easier fabrication process was intro- 
duced by J. Li et al in “Chemical and Physical Sensors, 
Carbon Nanotubes: Science and Applications”, M. Meyyap- 
pan, ed., CRC Press, Boca Raton, Fla., 2004. First, an inter- 
digitated electrode (“IDE”) configuration is fabricated using 
15 conventional photolithographic methods with a nominal fin- 
ger width of 1 0 pm and gap size of 8 pm. The electrode fingers 
are made of thermally evaporated Ti and Au (20 mn and 40 nm 
thickness, respectively) on a layer of Si0 2 , thermally grown 
on top of a silicon wafer. Second, a thin layer of carbon 
20 nanotubes forming a network is laid on the fingers using a 
solution casting process. The conductivity of the CNT net- 
work changes upon exposure of the fingers to different gases 
or vapors. Such a process is significantly simpler and pro- 
duces consistent sensors based on statistical properties of the 
25 CNT network with high yield (-100 percent). The IDE con- 
figuration facilitates effective electric contact between 
SWCNTs and the electrodes over large areas, while providing 
good accessibility for analytes in the form of gas/ vapor 
adsorption or contaminants adsorption/extraction in/from 
30 liquid to all SWCNTs including semiconducting tubes. 

Although carbon nanotube-based chemical sensors using 
the IDE configuration have been prototyped in laboratories, 
use of an IDE configuration for wireless sensing has not been 
reported in published research. Use of an IDE configuration 
35 requires frequent transmission of data measured at each of a 
large number of IDE fingers. A straightforward approach may 
consume substantial power and may require use of a relatively 
large footprint for the sensor system. This may be inconsistent 
with automated use of a remote system in a confined space. 
40 What is needed is a wireless transmission system that (1) 
has a small footprint and whose size is compatible with a 
nanosensor system, which may have a diameter as small as 
5-15 cm, (2) consumes relatively little power (e.g., about 50 
jli Watt-60 mWatt, or smaller with voltage regulation not acti- 
45 vated), (3) is reliable, with a mean time to failure exceeding 8 
hours (representative battery life) and (4) permits frequent 
data transmission from each of a laige number of data 
sources. 

50 SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides a 
system with one or more sensors for selected chemicals (all 
using nanostructure sensors with small physical sizes), one or 
55 more sensors for physical parameters, a multiplexer to receive 
and interleave the measured data stream values from the 
sensors according to a selected interleaving pattern, and a 
wireless transmission module to transmit the measured val- 
ues to a receiver and data analyzer. 

60 The overall sensor system consists of a chemical sensor 
module, a microcontroller-based data acquisition module, a 
multiplexer and constant current source module, and a wire- 
less communication module. The chemical sensor module is 
based on use of an interdigitated electrode (“IDE”) configu- 
65 ration. A layer of single-walled carbon nanotubes 
(“SWCNTs”) is laid on the IDE fingers. Chemical sensing is 
based on detection of changes in the conductivity of the 



US 7,968,054 B1 


3 

SWCNT network. In one embodiment, the system has 32 
channels of chemical sensing elements. It is self-contained 
and portable, and wirelessly transmits measurement data to a 
PC, using an I.E.E.E. 802.11a, 802.11b or 802.15 wireless 
LAN protocol. The footprint of the invention has a diameter 
as small as a few cm. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 schematically illustrates a sensor array and a wire- 
less transmitter according to an embodiment of the invention. 

FIGS. 2 A and 2B graphically illustrate sensor responses 
for identical exposures. 

FIG. 3 shows the variance of measurement data from the 
sensor in FIGS. 2 A and 2B, computed with a moving window 
of 1000 samples of the data. 

FIGS. 4 and 5 illustrate correlation of two sensors within a 
moving window. 

FIGS. 6 and 7 graphically illustrate relative chemical sen- 
sor response for exposure to N0 2 and for Cl 2 , respectively. 

FIG. 8 graphically illustrates separation of the gases N0 2 , 
Cl 2 , benzene and acetone according to principal component 
analysis. 

FIGS. 9,10 and 11 are flow charts of algorithm procedures 
for estimating concentrations of initially unknown constitu- 
ents according to the invention. 

FIG. 12 graphically illustrates transition from linear to 
logarithmic dependence upon a variable, such as concentra- 
tion. 

DESCRIPTION OF BEST MODE OF THE 
INVENTION 

FIG. 1 schematically illustrates components of a system 
10, constructed according to the invention, including a chemi- 
cal and/or physical nanosensor module 11 , whose output 
signals are received by a multiplexer (“MUX”) 12 . An output 
signal of the MUX 12 is augmented by a constant current 
source 13 and is passed through a buffer amplifier 14 (op- 
tional) and is received and processed by a microcontroller 15 
(e.g., a PIC16F688, available from Microchip Technologies). 
In one embodiment, wherein an array of carbon nanostruc- 
tures is applied to environmental sensing, the microcontroller 
15 also receives a signal from at least one of a temperature 
microsensor 16 , a vapor pressure microsensor 17 , a flow 
microsensor (direction, velocity) 18 , a relative humidity 
microsensor 19 and an acoustic or electromagnetic distance 
sensor 20 . The microcontroller output signals (e.g., inRS-232 
format) are received by a wireless communication module 21 , 
which includes a wireless server 22 (e.g., a MOXA W2150), 
which transmits these signals, using 802.11b protocol, 
through an antenna 23 to a local area network (LAN). The 
transmitted signals are received by a receiver (not shown) that 
is spaced apart from the transmission system 11 , for display 
and/or further signal processing. 

Each of the individual sensors, such as 16 , 17 , 18 and 19 , 
has its own data reporting cycle, and it is assumed herein that 
these reporting cycles are numerically compatible. In a first 
approach, each of the reporting cycles has the same length At, 
and the four sensors report to a multiplexer in a consecutive 
interleave pattern as t m fJ ={16, 17, 18, 19, 16, 17, 18, 19, 16, 
17, ... } (n=l,2, 3, 4;m=l, 2, . . .) In a second version, at least 
one of the sensors ( 16 ) has a reporting cycle length that is N 
times as long as the other cycle lengths (e.g., N=3), and the 
four sensors report in a second interleave pattern as t m „={16, 
17, 18, 19, 17, 18, 19, 17, 18, 19, 16, 17, 18, 19, 17, 18, . . . }. 
In a third version, the reporting cycle lengths are rational, 
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non-integer, multiples of each other so that the reporting 
pattern appears random, for example, where At(16)=(%)At 
(17)=( 3 A)At(18)=( 7 /4)At(19), a consistent third interleave pat- 
tern is: t mn ={\6, 18, 16, 17, 18, 18, 19, 16, 18, 18, 16, 17, 18, 
5 19, . . . }. For any interleave pattern w for data report- 

ing, the multiplexer transmits the received data in another 
time sequence {t 1 } , where f is determined with ref- 

erence to t m ^ (e.g., fl m w =t m)W +At, where At is a substantially 
constant time delay value). 

to In one embodiment, the chemical sensor module contains 
32 sensing elements, arranged in an IDE configuration, hav- 
ing pure SWCNTs, having polymer coated SWCNTs, and/or 
having metal nanoclusters or doped SWCNTs. At the center 
of the data acquisition system is a microcontroller 
15 (PIC16F688) from Microchip Technologies Inc. that samples 
each sensor element through a set of four Texas Instruments 
CD4051 multiplexers. Each MUX 12 reads signals from a 
group of eight chemical sensing elements. The LM234 con- 
stant current source from National Semiconductor is used to 
20 provide a constant current (100 pAmp) to each sensing ele- 
ment. Four of these devices are used to excite each group of 
eight chemical sensing elements. Conductivity or resistance 
is measured by supplying a constant current and measuring 
the corresponding voltage difference across the sensor. Also 
25 included in the data system is an AD22100K temperature 
sensor from Analog Devices. The microcontroller reads all 32 
chemical sensor and temperature values and generates a serial 
data output that can be connected directly (1) to a wireless 
serial device server, for wireless data transmission, or (2) to 
30 an RS-232 serial data output for a PC, for data logging. In this 
design, a W2150 serial server from Moxa Technologies, Inc. 
is used to wirelessly transmit data to a laptop or to a desktop 
PC, using an I.E.E.E. 802.1 lb wireless LAN protocol. 

An alternative server is a Crossbow MICAz, using an 
35 I.E.E.E. 802.15.4 protocol at 2.4 GHz, can also be used for 
wireless transmission. This system has a data transport rate of 
250 kbits/sec, and each node can also serve as a router. This 
allows a subset of one or more nodes to serve as a collection 
point for collecting data that were measured elsewhere and 
40 routing the data to a data processing module and/or to a 
control site that can respond to any problem encountered in 
the data measurements or data processing. This system uses 
direct sequence spread spectrum radio transmission, which is 
useful where two or more data collection nodes are transmit - 
45 ting substantially simultaneously to a data processing module 
and/or to a control site. Transmission range is up to 100 M 
outdoors and up to 30 M indoors. Device dimensions are 
58x32x7 mm (diameter «68 mm). Analog inputs as well as 
digital inputs are accepted by the MICAz system. 

50 Another alternative server is an Air-Pro D-Link dual band 
wireless LAN, operating at 5 GHz and 2.4 GHz, using 
801 .1 la/b protocol with 64/128 and 152-bit, shared key data 
encryption for authentication and enhanced data security. 
Each node can also provide data collection and forwarding to 
55 a data processing module and/or to a control site. Data trans- 
port rates are 1,2, 5.5, 1 1 Mbits/sec (802. lib) and 6, 9, 12, 18, 
24, 36, 48 and 54 Mbits/sec (802. 1 la). Transmission range is 
up to 400 M outdoors and up to 100 M indoors. Device 
dimensions are 118x54x13 mm (diameter «T33 mm). 

60 FIG. 1 schematically illustrates components of the system 
10, which includes a chemical nanosensor module 11 , whose 
sensor output signals are received by the MUX 12 . The output 
signal from the MUX 12 is augmented by a constant current 
source 13, and this signal is passed through a buffer amplifier 
65 14 (optional) and received by a microcontroller 15 (e.g., a 
PIC16F688, available from Microchip Technologies). The 
microcontroller 15 also receives a physical sensor signal from 
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at least one of a temperature micro sensor 16, a vapor pressure 
microsensor 17, a flow microsensor (direction, velocity) 18, a 
relative humidity microsensor 19 and an acoustic or electro- 
magnetic distance sensor 20. The microcontroller output sig- 
nals (presented, e.g., in RS-232 format) are received by a 5 
wireless server 21 (e.g., a MOXA W2150), which transmits 
these signals, using an I.E.E.E. 802.1 lb protocol, through an 
antenna 22 and are received by a receiver and analyzer 30, 
spaced apart from the transmission system 10 , for display 
and/or further signal processing. 10 

Sensor Response 

FIG. 2A shows typical sensor responses with the calcula- 
tion. FIG. 2B shows the sensor responses from three identical 
sensors that were observed over a period of 25 hours from 
three sensing elements . The graph i llustrates change of sensor 1 5 
resistance changes in response to the introduction of an ana- 
lyte, in this case ammonia, at approximately t=19 hours and 
t=20.5 hours. Data from the sensors were collected at one 
sample per second. For this experiment, a cotton swab soaked 
in ammonia hydroxide solution was placed approximately 2 20 
cm above the sensor surface, but the concentration of the 
analyte was not quantified. The sensor resistance for sensors 
x, y, and z increases approximately 1 0 percent from the sensor 
resistance measured just prior to the time(s) of detection of 
the analyte. Statistical properties of the measurement data 25 
were investigated to extract the information that presents the 
detection of analytes. Such information can be drawn from 
measurement data of a single sensor or from multiple sensors. 

FIG. 3 shows the variance of the measurement data from 
sensor “x” computed over a moving window of 1 000 samples 30 
of the data. This window represents 1 6.7 minutes of measure- 
ment data. Small variance indicates small changes in mea- 
sured data values over this period of time. A large variance 
indicates possible presence of analytes. In FIG. 3, this vari- 
ance is substantially larger at two instants of time (t=19 and 35 
20.5 hours) when the analyte is introduced. 

In addition to computing variance of individual sensor data 
values, statistical correlation of two sensors can be computed 
to confirm the presence of analytes. FIG. 4 shows the corre- 
lation coefficient for two sensors y and z, computed over the 40 
same moving window. FIG. 4 indicates, from an agreement 
between two sensors, that both sensors were exposed to the 
presence of analytes. An enlarged view of FIG. 4 around the 
two times of analyte introduction is shown in FIG. 5. Similar 
results were also observed for sensor pairs (x,y) and (x,z). 45 
FIGS. 6 and 7 graphically illustrate relative sensor response 
for exposure of the chemical sensor array to N0 2 gas and to 
Cl 2 gas, respectively. The functionalized sensor-to- sensor 
response differences can be used to determine or estimate the 
gas that is present. 50 

Chemical Discrimination 

The nanomaterials used in this array were individually 
studied for different gases and can be classified as: 1) pristine 
single walled carbon nanotubes, 2) carbon nanotubes coated 
with different polymers (here, chlorosulfonated polyethylene 55 
and hydroxypropyl cellulose), and 3) carbon nanotubes 
loaded with Pd nanoparticles and monolayer protected clus- 
ters of gold nanoparticles. 

When exposed to a vapor-phase analyte, each sensing ele- 
ment in the array responds uniquely. A reproducible combi- 60 
nation of resistances, or “smellprint”, for each vapor/gas is 
manifest. The sensor response is measured as a bulk relative 
resistance change (alternatively, current change or voltage 
change): AR/R 0 , as in FIG. 3 A. To calculate the relative 
sensor response (SR), the initial sensor resistance R 0 (or cur- 65 
rent/ voltage), and the resistance R (or current/voltage) during 
the exposure are measured, and a normalized change, 


SR = 


Rt-Ro 

Ro 


( 1 ) 


is computed. 

By comparing the results of several gases, it is possible to 
determine which sensors will respond in which manner (i.e. 
large or small amplitude, positive or negative). Looking at the 
histogram of sensor responses from 32 sensors to different 
gases/vapors, patterns for each gas/vapor can be established 
for discrimination between tested gases/vapors. When the 
sensor responses of an unknown subject gas are compared 
with the smellprints of several known substances, the 
unknown subject gas can be identified by matching its pattern 
to, or minimizing a pattern difference relative to, one of the 
known substances in the library. 

However, the sensor responses often vary with concentra- 
tion of the gas, temperature, and other external factors. To 
account for this, the data gathered by the sensors is normal- 
ized so that the relative responses can be compared, in order 
to make the sensor information more accurate. In normaliza- 
tion, the length of all data vectors becomes the same by 
making the sum of the squares of vector components equal to 
a constant so that the vectors become fixed length vectors. 
The data were normalized using: 


£4 


= C\ . 


( 2 ) 


Here, k designates the sensor, i identifies the gas, and NV is 
the total number of sensors (32 in one embodiment). Theo- 
retically, c z can be any positive number, such as 1 . 

In addition to analyzing the relative responses, the data are 
autoscaled to unit variance and zero mean, by mean-centering 
followed by dividing by the standard deviation: 


4 = 


X\k Xfc 

Sk 


( 3 ) 


where x\ k is the autoscaled response, x ik is the relative sensor 
response, x k is the mean value of normalized response for that 
specific sensor, and s* is the standard deviation: 


s k = 


ivpriT 


(■ Xik - Xk) 


( 4 ) 


where NP is the number of independent responses. Auto scal- 
ing removes any inadvertent weighting that arises due to use 
of arbitrary physical units. In the absence of any information 
that would preclude its use, use of auto scaling is preferable 
for most applications. 

The processed data set contains the sensor responses from 
all 32 channels. Pattern recognition algorithms are very pow- 
erful tools to deal with a large set of data. For example, 
principal component analysis (PCA) was applied to the data 
set that was collected from this carbon nanotube based sensor 
array. The PCA is to express the main information in the 
variables X={x Ar lk=l, 2 . . . K)} using a smaller number of 
variables 

T=(t i ,t 2 ...t A }{A<K), ( 5 ) 
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referred to as principal components of the variable X. In this 
instance, the X is a matrix of 32 sensor responses to these 
gases and vapors from the sensor array. Here, three principal 
components are labeled as FI, F2, and F3 that correspond to 
1 5 2? andt 3 in Eq. (5). FIG. 8 illustrates separation of the gases 5 
N0 2 , Cl 2 , benzene and acetone in terms of these three prin- 
cipal components. 

These gases and vapors are completely separated in prin- 
cipal component space, which indicates that this sensor array 
provides a high discrimination power to these gases and 10 
vapors. Because these gases were tested by this sensor array 
at different concentration levels in the range of 5-45 ppm, this 
sensor array can be used (uniquely) for high sensitive gas and 
vapor detection and discrimination. This sensor array dis- 
criminates the gases and vapors by their chemical nature 15 
rather than their concentrations. Also, the array responses can 
discriminate between gases that have some similarity in their 
chemical nature. 

Further analysis is often required to reliably estimate con- 
centration of one or more target chemicals cO. After exposure 20 
of a nanostructure to a gas, a measured electrical parameter 
value EPV (e.g., impedance, conductivity, capacitance, 
inductance, etc.) changes with time in a predictable manner, if 
a selected chemical precursor is present, and will approach an 
asymptotic value promptly after exposure to the target chemi- 25 
cal. The measured EPVs are compared with one or more 
sequences of reference EPVs for one or more known target 
precursor molecules, and a most probable concentration 
value is estimated for each of one, two or more target mol- 
ecules. An error value is computed, based on differences for 30 
the measured and reference EPVs using the most probable 
concentration values. Where the error value is less than an 
error value threshold, the system concludes that the target 
molecule is likely. Presence of one, two or more target mol- 
ecules in the gas can be sensed from a single set of measure- 35 
ments. 

For relatively large concentrations, it is assumed that the 
measured EPV for a given target molecule will vary linearly 
with the concentration C, 

40 

EP V(C)=a+b-C, tm ( 6 ) 

where the parameters a and b are characteristic of the particu- 
lar target molecule for which the measurements are made. 
Where EPV(C M ) and EPV(C 12 ) are values for a pure sub- 
stance measured at the different concentration values C l ?1 and 45 
C l 2 , respectively, the parameter values a and b can be esti- 
mated by 

a={EPV{C ul )C ia -EPV{C i;2 )C Ui }/{C ia -C Ui }, (7-1) 
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which CNS sensor sub -arrays are more sensitive to presence 
of a particular target molecule. In a first embodiment, the 
array is exposed to a selected target molecule, such as H 2 0 2 , 
at a selected sequence {C } (q=l, . . . , Q with Q=N) known 
(not necessarily distinct) concentration values (e.g., 500 ppm, 
14,000 ppm, 65 ppm, 1200 ppm, 10 ppm, etc.), and a refer- 
ence measurement value EPV(n;q;k;ref) (n=l, 2, . . . , N) is 
recorded for each sensor sub-array n, each concentration C 
of a reference gas containing a selected target molecule (k). 
Measurements for one or more concentration values for the 
selected target molecule gas can be repeated, if desired, to 
obtain an NxN square matrix of values. With q fixed, an NxN 
matrix {EPV(n;q;k;ref)/EPV(n;k;norm)}=E(n;q;k;ref) is 
formed, where, EPV(n;k;norm) is a normalization factor for 
the selected target molecule gas no. k and the sub -array no. n, 
which may be chosen as 

(n E/p (10) 

EPV{n\ k\ norm ) = < ^ u q EPV{n\ q ; k\ ref) p i , 

where {u^}^ is a selected set of non-negative weight numbers 
whose sum is a selected positive number (e.g., 1 or N) and p 
is a selected positive number. The normalization factor E(n; 
k;norm) may be a single EPV (e.g., u^=l andu^=0 for q^ql), 
or may be a weighted sum of two or more EPVs. This nor- 
malization (optional) is intended to compensate for the con- 
centration dependence of the particular target molecule. 

For each sub-array n (fixed), the mean and standard devia- 
tion for each of K reference gases, numbered k=l, . . . , K 
(Ki^2) are computed, as 

N (11 A) 

ju(«; k) = ^ E{n\ q\ k\ ref)/ N, 

<7= 1 

k (1 IB) 

fi(n) = J^fi(n; k)/K, 

*=1 

[n I 1 / 2 (12A) 

cr(«; k) = 1 2] {£(«; q\ k\ ref ) - / u(n- k)} 2 / A j , 

n (12B) 

cr(re) = ^ cr(n; k) / K. 

q=l 


b=(EPV(C ul )-EPV(C ia }/{C ia -C ul }, (7-2) 

for each target molecule. For smaller concentrations («l-50 
ppm), it may be preferable to use a logarithmic approxima- 
tion, 

EPV(C)=a'+b'-log e C, (8) 

where a' and b' are selected parameters. The parameter values 
a' and b' can be determined in a manner similar to that of Eqs. 
(7-1) and (13-2), by replacing the variables C M and C l 2 by 
log e {C 1 sl } and log e {C 12 }, respectively. FIG. 12 graphically 
illustrates the two concentration regimes for EPV(C), corre- 
sponding to Eqs. ( 6 ) and ( 8 ), which join together at a concen- 
tration transition value C tr for which 

a'+b'-io^C^a+b-C^. (9) 

Given an array of N sub-arrays of CNS sensors and a set of 
K target molecules (k=l, . . . , K;K^2), the EPV data for the 
sensors can be pre-processed in order to identify more clearly 


One now forms K autoscaled NxN matrices, defined by 

S(n;q;k)={E(n;q;k;ref)-[i(n)}/o(n), (13-k) 

and analyzes K eigenvalue equations 

S{n;q;k) V(k;Mk))=Mk) V(k;K(k)), (14-k) 

where k (=1, . . . , K) is fixed and V(k;X(k)) is a normalized 
Nxl vector that will usually depend upon the reference gas 
(k). If, as is likely, the N eigenvalues X(k) for a fixed reference 
gas k are distinct, the corresponding eigenvectors V(k;X,(k)) 
are mutually orthogonal (non-degeneracy). In the unusual 
event (degeneracy) that two or more of the N eigenvalues 
X(k)are equal, different non-zero linear combinations of the 
corresponding eigenvectors V(k;X(k)) can be constructed that 
are orthogonal to each other, within a sub -space spanned by 
the reduced set of eigenvectors corresponding to the identical 
eigenvalues. 


50 


55 


60 
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Each matrix equation (14-k) has a sequence of N (eigen- 
value;eigenvector) pairs, {(^(k);V M (k;X (k))}„, for a fixed 
reference gas k, and it is assumed here that the eigenvalues are 
arranged so that 

l*. 1 (*)l^l\ 2 (*)^ ■ ■ ■ =1^)1 (15-k) 

and so that the highest magnitude eigenvalue in each set 
satisfies 

IX 1 (*l)I^IX 1 (*2)l^ . . . IXi(iWV)U (16) 

where {kl, k2, . . . , kN} includes each of the integers {1, 
2, . . . , N} precisely once. The eigenvector V(kl;X 1 (kl)) is 
identified as a first basis vector V'(l): 

P(kl;X 1 (A:l))=r(A:l;X 1 (A:l)) (17-1) 

A second modified vector 

V'(k2 ^{k2))=V{k2 Ai(k2))-{ V(k2;X 1 (k2)), V\k\ 

(kl))}V'(kl-Mkl)) 17-2) 

is computed, where {V(k2;X 2 (k2)), V'(l)} is the scalar prod- 
uct (also referred to as the inner product) of the vectors 
V(k2;E 2 (k2) ) and V'(kl). More generally, a pth modified vec- 
tor 


V'(kp-x l (kp))= (17- p) 

p - 1 

V(kp- X ,i(*p)) - y {V(kp- X i(kp)), V'(kr, Xi{kr)))V’ (kr, *,(*/•)) 


is computed for p=2, . . . , K. The set of vectors {V'(kr;X 1 (kr)) 
} t (r=l, . . . , K) is mutually orthogonal, in the sense that the 
scalar products satisfy 


{ V ' (, kr ; Ai (kr)), V' (i ks ; X l (ks))} = 6 r , s . (18) 

> 0(r = 5 ) 

= 0 (r ± s). 

Each of the set of vectors {VXkr;>q(kr))} r (r=l, . . . , K) is 
maximally independent of each of the other vectors in the set, 
in the sense of mutual orthonormality (Eq. (18)). Each vector 
V , (kr;X 1 (kr)) will have relatively large (primary) contribu- 
tions from some of the sensor sub-arrays and will have 
smaller (secondary) contributions from the remainder of the 
N sub-arrays. The vectors V , (kr;X 1 (kr)) identify a maximally 
independent set of linear combinations of EPV responses 
from the N sub -arrays that can be used to distinguish presence 
of one reference gas (target molecule kr) from presence of 
another reference gas (target molecule ks). For example, if the 
set of reference gases are H 2 0 2 , H 2 0 and CH 3 OH, N=3 and 
three matrix eigenvalue equations are to be solved in Eqs. 
(9-k) (k=l , 2, 3). More generally, presence or absence of any 
of K target molecules (Ki^2) may be estimated. 

The linear combinations LC(kp) of EPV measurements for 
the different sensor sub-arrays correspond to modified prin- 
cipal components for the particular reference gases chosen. 
Choice of another set of another set of reference gases will 
result in a different set of modified principal components, 
although change of one or more concentration values within 
a reference gas may have little or no effect on the modified 
principal components. 

The preceding analysis concerning EPV measurements 
and estimates of concentration can be done wholly at the data 
processing site, partly at each of the data processing site and 
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at the chemical sensor module site, or wholly at the chemical 
sensor module site. The system provides an integrated chemi- 
cal sensor system based on use of SWCNTs as sensing ele- 
ments, including some preliminary testing results. The asso- 
5 dated data transmission system is wireless, small, portable, 
and readily deployable in the field. 

A nanostructure, or assembly of such structures, can be 
grown, for example, by a procedure discussed in connection 
with FIG. 1 in “Controlled Patterning And Growth Of Single 
to Wall And Multi-wall Carbon Nanotubes,” issued to Delzeit 
and Meyyappan in U.S. Pat. No. 6,858,197, incorporated by 
reference herein. 

What is claimed is: 

1. A system for receiving, analyzing and communicating 
1 5 results of sensing chemical and/or physical parameter values, 

the system comprising: 

N sensors numbered n=l, . . . , N (Ni^2), of parameters, 
where each sensor senses an electrical parameter value 
EPV(t m „;n) of a selected chemical or a physical param- 
20 eter at a selected sequence of times {t^t m w } m (m=l, 
2, . . . ), where at least one sensor senses a data set that 
includes at least one of a local temperature, a local vapor 
pressure, a local relative humidity, at least one of an 
acoustic signal and an electromagnetic signal for esti- 
25 mating a distance between the sensor and a reference 
surface, and presence of a selected chemical in a test gas, 
and where each sensor comprises a nanostructure (NS), 
having a sensing element having a diameter no greater 
than about 20 nm; 

30 wherein at least one NS is loaded with a selected sensitiz- 
ing substance drawn from a group of sensitizing sub- 
stances comprising Au particles, located between first 
and second ends of the at least one NS, and the first and 
second ends are connected to first and second terminals, 
35 respectively, of at least one of a voltage source and a 
current source, and to an electrical parameter value mea- 
surement mechanism that measures a change in at least 
one of electrical current, voltage difference, electrical 
resistance, electrical conductance and capacitance 
40 between the first and second ends of the at least one NS ; 

a multiplexer, having at least N input terminals and at least 
one output terminal, that receives at least one of a sensed 
value EPV(t m „;n) and a sensed value change AEPV 
(t m w ;n) from sensor no. n, at a sequence of times, 
45 {t=t’ m w }, where the time t' m „ is determined with refer- 
ence to a corresponding time X , and where each of a 
first sequence of times is interleaved with a 

second sequence of times {t m „ 2 } m for at least two inte- 
gers, nl andn2, satisfying l^nl<n2^N; 

50 a wireless transmission module, connected to the multi- 
plexer output terminal, to receive and transmit at least 
one of the sensed value EPV(t m „;n) and the sensed value 
change AEPV(t m „;n), received from the multiplexer, for 
at least one sensor no. n; and 

55 a computer that is programmed to receive, from the wire- 
less transmission module, and to compare at least one of 
the second value EPV(t m „;n) and the measured sensed 
value change AEPV(t m „) in an electrical parameter 
value with a corresponding reference value and to esti- 
60 mate and indicate, based upon the comparison, presence 
or absence of at least one chemical in a group of K target 
chemicals, numbered k=l, . . . , K (K^l), 

where the N sensors, the multiplexer and the wireless trans- 
mission module are contained in a volume having a 
65 diameter no greater than about 1 4 cm. 

2. The system of claim 1, wherein said multiplexer is 
located adjacent to at least one of said sensors. 
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3. The system of claim 1, wherein: 
at least one of said sensors is exposed to said test gas to be 
interrogated for presence or absence of one or more of 
said K target chemicals in said test gas (k=l, . . . ,K; 
K^l); 

measurement values EPV(n;test) are provided of a selected 
electrical parameter for said test gas for each said of 
sensor no. n=l, . . . , N (Ni^2), for each of said K target 
chemicals and for each of a selected sequence of con- 
centrations Ck (k=l, . . . , K; K^N) of the target mol- 
ecule no. k; and 

wherein said computer is further programmed: 
to provide at least N+l functional relationships 


K 

EPV{G ref ; meas\ n ) = a^ n + a k>n Ck(ref), 

k=l 

EPV{G re f \ meas\ n\C- 0) = a 0>n , 

relating a reference gas configuration, G re ^G({Ck(ref)}), in 
which target chemical no. k, is present in a reference gas with 
a known, non-negative concentration value Ck(ref) 
(k=l , . . . , K), to a test measurement value EPV (G^meas ;n) 
of the reference gas that is sensed by said sensor no. n, where 
a 0 and a Av? , are known calibration coefficients that are inde- 
pendent of concentration values present, and a 0 is an EPV 
value that would be measured by said sensor no. n in a test gas 
when all of the concentration values Ck are 0; 

to receive or provide a sequence of test measurement val- 
ues, V w (G;meas) and V 0 (G;meas), on said test gas in a 
gas configuration G=G({Ck}), where at least one of the 
concentration values Ck of a target chemical is 
unknown; 

to provide a numerical-valued error function €({Ck}) asso- 
ciated with the test measurement values V„(G;meas), 
where the error function is defined by 


2 s({Ck}) = z- EPV(G ; meets ; n) - EPV{G\ meas\ 0) - ^ a^Ckl 


where w„ are selected non-negative weight coefficients and p 
is a selected positive number; 

to provide at least K equations relating the calibration 
coefficients a 0 k and a w A: for k=f, . . . , K, by partially 
differentiating the error function e({Ck}) with respect to 
each of the concentration values Ck, and setting each of 
the partially differentiated expressions of €({Ck}) equal 
to 0; and 

to obtain solutions Ck for the at least K equations, for 
k=l, . . . , K, and to interpret these solutions as optimal 
solutions for the concentration values Ck. 

4. The method of claim 3, wherein said computer is further 
programmed to choose said positive number p to be p=2P, 
where P is a positive integer. 
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5. The system of claim 3, wherein said computer is further 
programmed: 

to receive a test measurement value for each of said N 
sensors for a new test gas, where said concentration 
5 value of at least one of said K target molecules in the new 
test gas is not yet known; 

to provide said error function €({Ck}), where said calibra- 
tion coefficients, a 0 Ar and a w Ar , are determined as in claim 
3 and said concentration values Ck (k=l, . . . , K) are 
to treated as initially unknown; 

to minimize said error function €({Ck}) with respect to a 
choice of value of at least one selected concentration 
value C(k=kl) (f^kf^K); and 

to interpret the at least one selected concentration value 
15 C(k=kf) as a most likely value of said concentration 
value for said target molecule no. kl in said new test gas. 

6. The system of claim 5, wherein said computer is further 
programmed: 

to minimize said error function €({Ck}) with respect to a 
20 choice of value of each of said selected concentration 
values Ck (f =k^K) for said new gas; 

to compute a value of said error function €({Ck})=€(min) 
using said selected concentration values for said new 
gas; and 

25 when the value €(min) is no greater than a selected error 
function threshold, e(thr), interpreting satisfaction of 
this condition as indicating that target gases are the 
primary components of said new gas. 

7. The system of claim 1, wherein a frequency of transmis- 
30 sion of a first reporting cycle for said first sequence changes 

monotonically relative to a frequency of transmission of a 
second reporting cycle for said second sequence in response 
to change of a length of said first reporting cycle relative to a 
length of said second reporting cycle. 

35 8. The system of claim 7, wherein said computer is further 

programmed: 

to minimize said error function €({Ck}) with respect to a 
choice of value of each of said selected concentration 
values Ck (1 ^k^K) for said new gas; 

40 to compute a value of said error function €({Ck})=€(min) 
using said selected concentration values for said new 
gas; and 

when the value €(min) is greater than a selected error 
function threshold, €(thr), interpreting satisfaction of 
45 this condition as indicating that at least one gas that is not 

one of said target gases is also a primary component of 
said new gas. 

9. The system of claim 1, wherein said frequency of trans- 
mission of said first reporting cycle decreases monotonically, 

50 relative to said frequency of transmission of said second 
reporting cycle, as said length of said first reporting cycle 
increases relative to said length of said second reporting 
cycle. 

10. The system of claim 1, wherein said sensors, said 
55 multiplexer and said wireless transmission module together 

consume no more than about 60 mWatt electrical power in 
operation. 



